www.mwe-journal.org 



Journal of Mining World Express Volume 2 Issue 4, October 2013 



Experimental Study of Compaction Creep 
Model of Broken Rock 

Shuncai Li* 1 ' 2 , Qiang Zhang 3 , Zhanqing Chen 2 

1 School of Mechanical and Electrical Engineering, Jiangsu Normal University, Xuzhou 221116, Jiangsu, China 

2 State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and 
Technology, Xuzhou 221008, Jiangsu, China 

3 Department of Aerospace Ground, Air Force Logistics College, Xuzhou 221000, Jiangsu, China 
n zsclsc@263.net; 2 215561271@qq.com; 3 chenzhanqing@vip. 163.com 



Abstract 

By MTS rock mechanics test system, the creep experiments at 
five stress levels were carried out for the water-saturated 
broken rock and the natural water-content broken rock with 
the same mass and different initial stacking density, and the 
displacement-time curves corresponding to these two 
samples have been obtained. Fitting the strain-time curves 
with Kelvin-Volgt creep model, we have got the creep 
parameters at each stress level. The results indicate that some 
factors such as compressive stress, initial stacking density 
and hydrous have distinct effect on the creep parameters. At 
lower stress level, owing to the existing of the internal pore 
water, the instantaneous elasticity modulus of the saturated 
broken rock is slightly lower than that of the natural 
water-content broken rock; while at the higher stress level, 
with the discharge of pore water, the instantaneous elasticity 
modulus of the sample with the bigger initial stacking 
density is also higher. Under the same stress level, the limit 
creep modulus of the natural water-content broken rock is 
higher than that of water saturated broken rock, but the 
viscosity coefficients of the these two samples have little 
difference. 
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Introduction 

In mining engineering, the compaction deformation 
properties of broken rock caused by mining have a 
great significance in predicting the surface subsidence 
and in analyzing the shifting law of the overlying 
strata and migration law of liquid and gas in the 
broken rock. The compaction deformation of the 
broken rock is affected by many factors such as the 
compressive stress, the rock category, the breakage 
degree, stacking density and hydrous state, which 
leads to the great difference in compaction property of 
broken rock. 

So far, a lot of research has been done on the 



compaction of the broken rock. Miao (1997) studied the 
swell and compaction properties of the broken coal 
and obtained the curves of compressive stress 
changing with swell coefficient and of the compaction 
deformation changing with lateral stress; and Zhang 
(2005) drew a conclusion by the compaction 
experimental study of the granular rock that the 
breaking rate is the biggest when the relative pressure 
is between 0.5 and 0.6, and the grading curve after 
crushing can be fitted with a quadratic polynomial. By 
the compaction experiment of broken coal, shale and 
sandstone, Ma (2005) obtained that the curve of the 
compressive stress and strain can be fitted by a 
exponential relationship; Jiang(2001) studied the 
compaction mechanism of broken coal gangue and 
pointed out that the compaction course can be divided 
into two stages: crushing compaction stage and 
consolidation compaction stage; while Guo(2002) 
studied the deformation property of cracked rock in 
goaf and obtained the distribution law of swell 
coefficient for the cracked rock in goaf during long 
wall mining process; and Xu (2011) studied the 
deformation properties of backfill body during 
gangue-fly ash filling mining and obtained the best 
ratio of gangue-fly ash backfill body. 

The above studies did not involve the time-dependent 
deformation of broken rock. In fact, the subsidence 
above the gob is a slow process which needs to deal 
with the rock creep. According to statistics, mining 
subsidence lands of the medium-sized state-owned 
enterprises in China get up to 84201.4 hm 2 [Guo, 2004]. 
During mining, with the periodic breaking and caving 
of overlying strata, the caving broken rock in goaf 
transit from the natural accumulation state to 
compaction state, and even if through a long time of 
natural compaction, the broken rock will experience 
the second deformation under all kinds of internal and 
external factors (such as ground stress, ground water 
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and ground construction loads, and etc). This process 
is called as old gob activation and can cause a serious 
threat to the ground buildings, so the gob management 
often encounters the compaction, seepage and creep 
problems of broken rock. In recent years, with the 
promotion and application of green mining technology, 
the studies about effects of these factors (such as 
compaction degree of gangue filling, the long term 
deformation and seepage) on the controlling ground 
subsidence have become the hot topics in mining 
engineering [Miao, 2007; Zhang, 2007; Qian, 2008]. 
According to the existing references, many 
experimental studies have been done on the creep 
properties of the intact rock [Zhu, 2002; Fan, 2005], but 
the research on the compaction creep of broken rock 
(especially with pore pressure or seepage) is very 
limited and it is mainly concentrated on deformation 
experiment study of dry rock samples [Cheng, 2004; 
Liang, 2003], for example, Shen (1998) proposed a 
three-parameter model of rockfill creep and Cheng 
(2004) gave a nine-parameter model. The slow 
deformation of broken rock can lead to the change in 
its porosity, and Ma(2009) studied the porosity 
changing with the time during the creep for the broken 
mudstone with different grain sizes, but failing to give 
the creep constitutive equation of mudstone. 

Based on the above studies, in this paper the 
compaction creep experiments under five stress levels 
for saturated broken rock and natural water content 
broken rock with grain size 10-15 mm were carried out 
by the MTS test system and compaction/seepage 
patented device, as a result, the stress-strain 
constitutive relations during the compaction creep of 
broken rock were established. 

Test Equipment and Rock Sample 

As shown in Fig. 1, the test equipment is MTS815.02 
rock mechanics test system. Rock sample is 10-15 mm 
broken sandstone with mass 1700 g. The uniaxial 
tensile strength of the sandstone core is 57.5 MPa, and 
its density P is 2473 kg/m3 at natural water content. 




Note: S1-S15 -switches; Ri~R2- reservoir; PT-pressure transducer; 
DPT-differential pressure transducer; I-pressure intensifier; 1-upper 
platen for axial loading; 2-piston; 3-cylinder tank; 4-felt; 
5-permeable plate with many holes; 6-rock sample; 7-base plate 

FIG. 1 MTS 815.02 ROCK MECHANICS TEST SYSTEM 
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Test Steps 

■ Determine the initial height 

Taking the rock sample with m=1700 g into the 
cylinder tank, according to the height H=180 mm, 
the total height of piston 110 mm, the height hi that 
the piston head outdoes the cylinder tank, the 
thickness of the permeable plate h2=9 mm and the 
thickness of felt h3 =3 mm, we can get the initial 
height ho 

h Q =180-(110-/z 1 )-2x(9 + 3) = 56 + /z 1 (mm) (1) 

■ Saturate the rock sample 

Before the axial loading, the rock sample must be 
saturated by the test system, but for the natural 
water content broken sandstone this step does not 
require. 

■ Axial force loading 

Exert each axial compressive force one by one, all 
the loading time of each force level is 40 s and the 
force is not unload when loaded to the next force 
level. After each level of force has been loaded, it is 
kept steady, and then rock samples started to creep. 
Each loading level and its keeping time are shown 
in Table 1. The MTS system was set up to collect the 
test data (including the displacement, loads) of rock 
sample at every 5 s. 



Table 1 Each Force Level And Its Keeping Time 



load/ kN 


Stress /MPa 


Keeping time of each force /min 


25 


cr, = 1.97 


15 


50 


cr 2 = 3.94 


30 


100 


<T 3 = 7.98 


60 


200 


<T 4 = 15.79 


120 


400 


a- = 31.58 


180 



■ Calculate the macro-strain (logarithmic strain) 

According to the collected axial displacement S, we 
can calculate the stacking height h and compressive 
strain s of sample at each moment t, 











£ = 






In(l-f-) 




K 




h 



(2) 



Fitting the s-t curve and a - s level, we can obtain 
the creep model parameters of broken rock. 

Test Results and Analysis 

Displacement-Time Curve of Saturated Broken Rock 

The mass m of saturated broken rock is 1700 g; whose 
initial stack height in the cylinder tank is 102 mm. 
Because the inner diameter D of cylinder tank is 127 
mm, we can calculate its initial stack density 
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4m 



= 1316 kg/m 3 . Displacement- time curves of 



Pp ~nD L 

saturated rock sample at 5 levels stress are shown in 
Fig. 2. 
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FIG. 2 DISPLACEMENT-TIME CURVES AT EACH LEVEL STRESS 

According to displacement-time curves, the deformation 
of broken rock tends to be stable after a primary creep. 



Creep curve has a distinct fluctuation at <J 1 =1.97 MPa, 

resulting from the fact that the position of sample 
grains suddenly changes under the axial compressive 
force. The deformations at each level stress include 
transient deformation at loading stage and creep 
deformation at stress keeping stage. 

Creep Model of Saturated Broken Rock 

The common creep models are as follows: Maxwell 
model, Kelvin model, Kelvin- Volgt model (K-V model), 
Burgers model and so on, as shown in Fig.3 (a)-(d). 
Creep in the Maxwell model and Burgers model 
increase with time and up to infinite. So Maxwell 
model and Burgers model are liquid models and 
inconsistent with actual situation of crushed stone, so 
they cannot be used. Kelvin model cannot reflect 
instantaneous elastic deformation of material, 
obviously, inconsistent with elastic properties of the 
test, so it cannot be used either. Only Kelvin- Volgt 
model can reflect creep law and instantaneous elastic 
property. Therefore, the creep of the broken rock has 
been analysed by Kelvin-Volgt model (K-V model) in 
this paper. 



^ ,/VV — ED °- 

(a) Maxwell model 
E 

A/V 



WW 1 



n 

3Zk 

(b) Kelvin model 
& 




(c) Kelvin-Volgt model 

s r , . , A/V — 

(d) Burgers model 
FIG. 3 TRADITIONAL CREEP COMBINATION MODEL 
The constitutive form of K-V model is 



o" cr t . _ 
*(*) = — + — [l-exp( )] 



(3) 



Where Eo is the instant elastic modulus, Ei is the limit 
creep deformation modulus, r i is delay time, and 



7i 



Where 1\ is viscosity coefficient. 



(4) 
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The instant elastic modulus Eo can be determined from 
stress and elastic strain at loading stage. The data 
collected in 30 s such as displacement S, axial force F at 
the first stress level a x = 1.975 MPa are as shown in 
Table 2. Fitting stress-strain diagram drawn by Table 2 
we can obtain Eo=13 MPa for the saturated broken 
sandstone at (7 X =1.97 MPa. 

Table 2 Stress And Strain At Loading Stage Of cr, 



t/s 


S/mm 


8 


F/N 


°"/MPa 

















5 


0.354 


0.00347 


2775 


0.219 


10 


0.895 


0.00882 


6061 


0.478 


15 


1.701 


0.01682 


8770 


0.692 


20 


3.591 


0.03584 


11673 


0.921 


25 


6.737 


0.06833 


15365 


1.213 


30 


9.386 


0.09653 


18831 


1.487 


35 


11.740 


0.1222 


21997 


1.737 


40 


13.712 


0.1443 


25272 


1.996 



Then we calculate creep parameters Ei 7i . Because 
instant elastic strain can be expressed as s = — , then 

the limit creep can be expressed as 



So we have 



— + — = + — 

E E x E x 



(5) 

According to the test data, for = 1 - 975 MPa, 
e = 0.1443 . When the time attains to 940 s, 
si s = 1.247 , Let £ «> = l-5£ 1 we can get Ei from Eq.(5) 

K = = — — = 27.35 MPa (6) 

1 0.5x0.1443 

From Eq.(3)we have 



s-s Q =-{\-e * ) 
E, 



(7) 



Then we have 



ln[cr - E 1 {s-£q )] = In cr -t 

Let y = \n[o--E ] (£-£ )] / a = -^- / b = \na, Then 
y = at + b 

By fitting test data, we have a = -0.1426 and 

//, =-^=191 MPa.s. 

a 



From the above, we have got Eo 12 MPa, E x 27 MPa 

and viscosity coefficient 1\ 191 MPa.s for saturated 
broken sandstone at the first level stress 
a x =1.975 MPa. Accordingly, we can obtain creep 
parameters of saturated broken sandstone at different 
stress levels as shown in Table 3, based on which, we 
have £ = 4.82a 1 1389 , E x = 10.05a 1 3029 , tj = 73.55a 1 349 . 

Table 3 Creep Parameters Of Saturated Sample 







a 2 


0- 3 


°4 


°" 5 


Eo /MPa 


13 


20 


40 


116 


281 


Ei /MPa 


27 


58 


135 


331 


1038 


711 /MPa.s 


191 


463 


1128 


3085 


7973 



The Displacement -Time Curve of Natural Water- 
Content Broken Sandstone 

The mass of natural water- content broken sandstone is 
m=1700 g, whose initial stacking height in cylinder 
tank is h =106 mm. The corresponding displacement- 
time curves (S-t curves) under five stress levels are 
shown in Fig. 4. 



45 

40 

35 
30 

e 

M 20 

15 

10 

5 




ff3=31J78MPa 

<j 4 =15.7S9MPa 

o. 3 =7.S95MPa 

o. 2 =3.947MPa 

CM=1.973MFa 



6iJi:nj 
tfs 



FIG. 4 S-t CURVES FOR NATURAL WATER-CONTENT SAMPLE 

The Creep Parameters of Natural Water-Content 
Broken Rock 

Using the calculation method in the previous section, 
we can obtain the creep parameters at different stress 
levels for natural water content broken rock as shown 
in Table 4, based on which we have 



E = 1.62a 



r ,E 1 = 23.97a 1197 , /7 = 213.23cr L 



Table 4 Creep Parameters Of Natural Water-Content Broken 
Sample 





<*i 


a 2 


0-3 


°4 


°" 5 


£0 /MPa 


19 


20 


37 


81 


200 


Ei /MPa 


50 


149 


269 


579 


1613 


/MPa.s 


460 


842 


1717 


3501 


8069 



According to Tables 3 and 4, the E -a curve, E x -<j 
curve and rj x -<j curve for the above two rock 
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samples are shown in Figs. 5-7. In Fig. 5, at the first 
stress level, the instantaneous elastic modulus Eo of the 
saturated broken rock is slightly lower than that of 
natural water-content broken rock, and the existing of 
pore water is the main reason leading to the decreasing 
Eo. With the compaction in progress, the pore water in 
saturated rock sample is drained out, the saturated 
rock sample is close to natural water-content broken 
rock; but because the initial stacking density of 
saturated sandstone is slightly bigger, at the same 
stress level, its compaction deformation is relatively 
smaller, so from the second stress level the 
instantaneous elastic modulus Eo of the saturated 
broken rock is slightly larger than that of natural 
water-content broken rock. From Figs. 6-7, it can be 
seen that with the increase of stress levels the limit 
creep modulus Ei and viscosity coefficient 7i also 
increase, and Ei of natural water-content broken rock is 
slightly bigger than that of saturated broken rock, but 
the coefficients 1\ of these two samples have little 
difference. 



300 r 




5 10 15 20 25 30 35 

a/MPa 

FIG. 5 E o-& CURVES OF TWO SAMPLES 
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FIG. 6 E l -a CURVES OF TWO SAMPLES 




5 10 15 20 25 30 35 

o/MPa 
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Conclusions 

By MTS rock mechanics test system, the creep 
experiments under five stress levels for 10-15 mm 
saturated broken rock and natural water-content 
broken rock were carried out, consequently, the 
displacement-time curves and Kelvin- Volgt creep 
model parameters were obtained. The research 
indicated that with the increase of the stress levels the 
creep parameters of the two samples increased 
markedly; at the same stress level, the saturated 
broken rock with a higher initial stacking density has a 
bigger instantaneous elasticity modulus and a lower 
limit creep modulus, but the viscosity coefficient of the 
these two samples has little difference. 
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